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Elastic and structural instability of cubic Sn 3 N 4 and C 3 N 4 under pressure 
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We use in-situ high pressure angle dispersive x-ray diffraction measurements to determine the 
equation of state of cubic tin nitride (7-Sn3N4) under pressure up to about 26 GPa. While we find 
no evidence for any structural phase transition, our estimate of the bulk modulus (-Bo) is 145 (± 
1.7) GPa, much lower than the earlier theoretical estimates and that of other group IV-nitrides. We 
corroborate and understand these results with complementary first-principles analysis of structural, 
elastic and vibrational properties of group IV-nitrides, and predict a structural transition of Sn3N4 at 
a higher pressure of 88 GPa compared to earlier predictions of 40 GPa. Our comparative analysis of 
cubic nitrides shows that bulk modulus of cubic C3N4 is the highest (379 GPa) while it is structurally 
unstable and should not exist at ambient conditions. 

PACS numbers: 61.05.C-, 62.50.-p, 62.20.-x, 63.20.D- 
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Since the discovery of cubic spinel (7) phases in Group 
IV nitrides [D-Q , they have received a renewed interest in 
the past few years due to prediction of super hardness [4j 
in these class of materials. It is also expected that these 
spinel classes of nitrides will exhibit interesting electronic 
properties, such as varying electronic band gap energies 
suitable for optoelectronic applications. So far, the bi- 
nary nitrides of tin[l[, germanium Q and silicon Q have 
been synthesized and shown to adopt the spinel structure 
under different experimental conditions. These new poly- 
morphs have the dense spinel structure with Si/Ge/Sn 
atoms in octahedral as well as tetrahedral coordination 
with nitrogen and the coordination of nitrogen is four. 
The crystal structure of 7-A3N4 (J-A3N4, A = Si, Ge and 
Sn) (space group Fd3m, #227) can be described as a dis- 
torted close-packed lattice of nitrogen atoms at (u,u,u) 
positions with A atoms occupying 1/8 of the tetrahe- 
dral interstitial sites and 1/2 of the octahedral sites [l]-[3j]. 
The anion arrangement produces a rigid-vertex linkage 
of A-filled regular tetrahedra (^4*N 4 ) and distorted oc- 
tahedral (^4°N6) sharing half of their edges. The cation 
sublattice can be described as an array of empty corner- 
sharing tetrahedra and A-filled truncated tetrahedral. 
The bulk modulus of the spinel phase is expected to be 
high, primarily determined by the nitrogen sublattice, 
but it should be influenced by the choice of cations. A 
very high bulk modulus of 369 GPaQ has been predicted 
for the cubic phase of C3N4, without any success in syn- 
thesizing this material. Though both Si3N 4 and Ge3N4 
have similar bulk moduli, the same for Sn3N 4 has been 
predicted to be significantly different. Also, from recent 
nanoindentation studies Q, tin nitride spinel was found 
to be significantly softer and more plastic than its lighter 
congeners (Si3N 4 and Ge 3 N 4 ). Thus, the exact role of the 
cation in determining the stability and elastic properties 



is not yet completely understood. 

Tin nitrides are of special interest due to their promis- 
ing semiconducting and electrochromic properties [71 |8j. 
The feasibility of utilizing tin nitride thin film as write- 
once, optical recording media has also been reported^. 
Therefore, the precise knowledge of the stability and me- 
chanical properties of 7-Sn3N4 would be of great interest 
for several practical applications. Though the ambient 
bulk modulus (Bq) of 7-Sn3N 4 and 7-Ge3N4 has been 



intensively investigated 1(3- l2] , there have been no ex- 
perimental studies to determine the compressibility of tin 
nitride. With lack of experimental data, first-principles 
calculations have been used to investigate the proper- 
ties of Sn 3 N 4 @,[I!-[l6|. 

However, the predicted values of 
the bulk modulus appear to vary between 186 and 218 
GPa. In this work, we have performed high pressure x- 
ray diffraction measurements on 7-Sn3N4 to determine 
the bulk modulus from the volume as a function of com- 
pression. We have also carried out first-principles den- 
sity functional theory (DFT) calculations in the gener- 
alised gradient approximation (GGA) to achieve a consis- 
tent understanding of our experimental findings. First- 
principles calculations have been carried out for other 
group-IV nitrides as well to (a) determine the stability 
of the spinel phase with respect to pressure and possible 
phase transitions, and (b) understand the influence of 
cation on elastic behaviour. Based on our calculations, 
we also discuss the stability of cubic phase of C3N4 at 
ambient conditions. 

Spinel tin nitride (Sn3N4) was synthesised by high- 
pressure solid-state metathesis reactions where tin 
tctraiodide was reacted with lithium nitride and ammo- 
nium chloride in a piston-cylinder apparatus at 623 K 
and 2.5 GPa. Details of the synthesis procedure can 
be found elsewhere 17j . In-situ high-pressure angle dis- 
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persive x-ray powder diffraction measurements (up to 
26 GPa) were performed using the monochromatic syn- 
chrotron radiation (A = 0.68881 A) at XRD1 beam- 
line at ELETTRA synchrotron source, Italy. The two- 
dimensional x-ray diffraction (XRD) patterns were col- 
lected on a MAR345 imaging plate. The sample-to- 
detector (image plate) distance was calibrated by col- 
lecting the diffraction pattern of powdered silicon (Si) 
at ambient conditions. The one-dimensional diffraction 
patterns were obtained by integrating along the Debye- 
Scherrer ring in the two dimensional image patterns using 
the FIT2D software [18j. Polycrystalline Sn 3 N 4 powder 
was placed together with gold into the 200 /xm hole of 
a stainless steel gasket (T301), preindented to 60 /im, 
inserted between the diamonds of a Mao-Bell-type di- 
amond anvil cell (DAC). Gold was used as an inter- 
nal pressure calibrant and the pressure in the DAC was 
determined using the known equation of state of gold 
[HSj . Methanol-ethanol-watcr (MEW) mixture (16:3:1) 
was used as pressure transmitting medium. To develop 
a microscopic understanding of the discrepancy between 
the calculated bulk modulus earlier [g, hJ-UJl and our ex- 
perimental estimate, we present new first-principles cal- 
culations based on DFT with a generalized gradient ap- 
proximation (GGA)[i3| for the exchange correlation en- 
ergy, as implemented in PWSCF[21| package. We use 
ultrasoft pseudopotentials[22| to describe the interaction 
between the ionic core and valence electrons, and a plane 
wave basis with energy cutoffs of 30 Ry and 240 Ry in 
representation of wave functions and charge density re- 
spectively. Integrations over Brillouin zone are sampled 
with an 8 x 8 x 8 Monkhorst Pack Mesh[23j. Struc- 
tural optimization is carried through minimization of 
energy using Hellman-Feynman forces in the Broyden- 
Flecher-Goldfarb-Shanno (BFGS) based method. Com- 
plete phonon dispersion is determined using DFT-linear 
response and 2x2x2 mesh-based Fourier interpolation 
of interatomic force constants. 

Figure 1 shows the XRD patterns of 7-Sn3N4 at var- 
ious pressures. The cell parameters were obtained by 
analysing the XRD profiles by Le-Bail profile fit using the 
GSAS software [iij]. Ambient lattice parameter has been 
found to be a = 9.0205(5) A, which is in close agreement 
to earlier reported value [ljj. Upon increasing the pres- 
sure, we didn't observe any new diffraction peaks (see 
Fig. 1) or sudden discontinuity in the pressure depen- 
dency of the d values up to the maximum achieved pres- 
sure of 26 GPa. This indicates that 7-Sn 3 N4 doesn't 
undergo any transition and remains in the cubic crystal 
structure up to the highest pressure. Figure 2 shows pres- 
sure volume compression data. Using the second order 
Birch-Murnaghan (B' = 4) equation of state (EOS) [25] 
to fit the volume compression data, we found the bulk 
modulus (B ) to be 145 (±1.7) GPa. It is interesting to 
note that the value thus obtained is much lower than the 
predicted bulk modulus values. Thus it was necessary to 
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FIG. 1. Pressure evolution of X-ray diffraction patterns of 
7-SnaN4 compressed quasi- hydrostatically in MEW up to 26 
GPa. Diffraction lines due to gold (Au) (used for in- situ 
pressure calibration) and stainless steel gasket (*) are also 
identified. 
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FIG. 2. Relative volume change in 7-Sn3N4 as a function 
of pressure at ambient temperature (Solid line is a second 
order Birch-Murnaghan fit to the experimental data (open 
symbols)). 



revisit the theoretical calculations. 

Our theoretical estimate of the lattice constant of 7- 
Sn3N4 is 9.136 A, within typical GGA errors of the exper- 
imental value of 9.0205 A[17j. Bulk modulus estimated 
from first-principles energies as a function of volume is 
158 GPa, in reasonable agreement with our measured 
value of 145 GPa re por ted here. We note that the earlier 
theoretical estimate [16J of the bulk modulus is within lo- 
cal density approximation (LDA) (contrary to the claim 
of use of GGA in this paper; we have reproduced their 
result with the same code and LDA choice of energy func- 
tional), and is overestimated to be 187 GPa. Thus, the 
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observed softness of S113N4 reflected in the lower value of 
B here should be reliable. With the same choice of ex- 
change correlation energy functional, we systematically 
studied properties of other group IV nitrides in the spinel 
structure, and find that the bulk modulus (see Table I) 
reduces from 379 GPa of carbon nitride to 158 GPa of 
Sn3N4. While the lattice constant increases from car- 
bon to tin by almost 30%, consistent with the size of the 
cation, the internal structural parameter (u) changes by 
only 2%. 

TABLE I. Experimental and calculated equilibrium lattice 
constant (ao) and bulk modulus (B) of the 7-phase of different 
group-IV nitrides. 

Calc.ao Expt.ao CalcB Expt.B u 

(A) (A) (GPa) (GPa) 

7-Sn 3 N 4 9.136 9.0144 158 145 (± 1.7) 0.2592 
7-Si 3 N 4 7.792 7.7381(2) 292 308(± 5) 0.2575 
7-C3N4 6.796 Not synthesized 379 0.2559 



The local stability of a structure can be confirmed 
through determination of the full phonon dispersion, and 
by verifying that there are no unstable modes in the 
structure. Our results (see Fig. 3) clearly show that 
C3N4 in the spinel structure is unstable with respect to 
both shear acoustic and optical modes, and should not 
form. This is consistent with the fact that there are no 
experimental reports so far on C3N4 in the spinel form. 
We also find that the frequency range (band-width) of 
acoustic branches reduces from silicon to tin nitrides, 
confirming their increasing elastic softness. If K is the 
spring constant or the stiffness of a bond and a is the 
lattice constant, a simple analysis shows that the bulk 
or elastic modulus should scale as K/a. Although the 
slopes of the acoustic modes in these nitrides seem dif- 
ferent in the phonon dispersion (see Fig. 3), variation of 
B with 1/a is almost linear. This means that the bond- 
stiffness in the nitrides studies here is quite similar, and 
it is the volume of the unit cell that is responsible largely 
for the monotonous decrease in the bulk modulus from 
carbon to tin nitride. As the spinel structure contains 
cations with tetrahedral and octahedral coordination, it 
is not straightforward to describe it with a single type of 
a bond. Since all anions are symmetry equivalent, our 
reference to bond stiffness here is reasonable in the sense 
of average bond stiffness. 

Our estimates of the Born effective charges also bear 
a similar feature (see Table II), with two types of effec- 
tive charges of a cation and degenerate effective charges 
of N. While the effective charges of SisN4 and Sn 3 N 4 are 
similar, those of C3N4 are singularly different. We find 
that the effective charge of carbon is most anomalous (i.e. 
deviates most from its nominal charge of Z*=4); this is 
not very surprising, as the anomalous charges (though 
larger than the nominal values) are known to correlate 
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FIG. 3. Phonon dispersion of a) 7-Sn 3 N4, b) 7-Si3N4 and c) 
7-C3N4 at equilibrium lattice constant. 



with structural instabilities (the presence of unstable vi- 
brational modes) in the system [26j. We also note that 
the effective charges of carbon atoms with 4-fold and 6- 
fold coordination are rather different, in contrast to rel- 
atively similar charges of the two symmetry inequivalent 
cations in Si and Sn nitrides. We find that the band 
gap steadily increases from 0.2 eV in Sn 3 N 4 to 3.4 eV in 
Si3N4, but reduces substantially to 1.02 eV in C3N4, con- 
firming rather different nature of the electronic structure 
and bonding in C3N4 from the rest. Our analysis based 
on projection of energy eigen-functions on the atomic or- 
bitals shows that carbon in spinel-C3N4 is close to being 
sp 3 -hybridized with covalent bonding with N, while the 



TABLE II. Calculated Born effective charges of the octahe- 
dral, tetrahedral coordinated cations and N atoms and elec- 
tronic dielectric constant in the 7-phase of different nitrides. 





Z* of Z* of 
tetrahedral cation octahedral cation 


Z* of N Dielectric 
constant 


7-Sn 3 N 4 


3.90 4.17 


-3.02 


9.17 


7-Si 3 N 4 


3.59 3.70 


-2.72 


5.29 


7-C3N4 


1.90 3.98 


-2.45 


8.99 
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FIG. 4. Phonon dispersion of 7-S113N4 as a function of pres- 
sure 



lus of 158 GPa, well below the earlier calculations. Thus, 
we demonstrate a significant increase in the compress- 
ibility for 7-Sn3N4, compared to lighter members (Si and 
Ge) in the family. We predict a structural phase transi- 
tion in 7-Sn 3 N 4 at 88 GPa based on phonon dispersion 
calculation as a function of pressure. While our calcula- 
tions confirm a very high value of bulk modulus of C3N4 
(as predicted earlier) [5j, our finding of unstable phonon 
modes and anomalous Born effective charge show that 
7-C3N4 is unstable at ambient conditions. This is in 
agreement with no experimental reports so far on C3N4 
synthesized in the cubic form. 

GKP, SKD and CN thank ICTP for the financial sup- 
port under the ICTP-Elettra users programme to carry 
out this experiment at Elettra, Trieste, Italy. GKP 
thanks Dr. Diptikanta Swain for useful discussions. SKD 
thanks Dr. M. P. Shemkunas for providing the samples. 



cations in Si 3 N4 and Sn 3 N4 appear more ionic. In sum- 
mary, the unexpectedly low measured bulk modulus of 
Sn 3 N4 does follow the trend in properties of the group 
IV-nitrides and seems to be readily understandable. 

Our experiments show that Sn 3 N4 remains stable in 
the spinel structure up to 26 GPa. We now use first- 
principles simulations to explore its stability at much 
higher pressures. Since it is difficult to explore all pos- 
sible structures as a function of pressure, we determine 
complete phonon dispersion of cubic Sn 3 N4 as a func- 
tion of pressure. This is an efficient method to find un- 
stable modes and detect any structural transition that 
the spinel structure would undergo with pressure, and is 
similar in spirit to approach used commonly in ferroelec- 
tric structural transitions [2 7j. Our results (see Fig. 4) 
show that the Sn 3 N4 remains stable locally in the spinel 
structure up to pressures of 88 GPa, and develops an in- 
stability in the transverse acoustic branch along <111> 
direction. Thus, Sn 3 N4 in the spinel structure will distort 
with a shear-strain above 88 GPa suggesting a transition 
from cubic to low symmetry phase. At higher pressures, 
we find that many more unstable modes develop that in- 
volve optical phonons as well, suggesting that the struc- 
tures beyond 88 GPa are expected to be more complex. 
Electronic structure of Sn 3 N4 is also rather sensitive to 
pressure, but exhibits an unusual trend: its band gap in- 
creases rather sharply from 0.2 eV at pressure to about 
2 eV at a pressure of 50 GPa. This is indeed consistent 
with the trend in band gaps from Sn 3 N4 to Si 3 N4, where 
the pressure has chemical origin. 

In summary, in-situ angle dispersive high pressure x- 
ray studies reveal that bulk modulus of 7-Sn 3 N4 is 145 
(±1.7) GPa and is stable up to a pressure of 26 GPa. 
Our first-principles calculations are in excellent agree- 
ment with these results with an estimate of a bulk modu- 
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